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ABSTRACT 

The Hubble Deep Field images have provided us with a unique chance to relate 
statistical properties of high-redshift galaxies to statistical properties of Lya absorption 
systems. Combining an empirical measure of the galaxy surface density versus redshift 
with an empirical measure of the gaseous extent of galaxies, we predict the number den- 
sity of Lya absorption systems that originate in extended gaseous envelopes of galaxies 
versus redshift. We show that at least 50% and as much as 100% of observed Lya 
absorption systems of W ^ 0.32 A can be explained by extended gaseous envelops of 
galaxies. Therefore, we conclude that known galaxies of known gaseous extent must 
produce a significant fraction and perhaps all of Lya absorption systems over a large 
redshift range. 
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1. INTRODUCTION 

Comparison of galaxy and absorber redshifts along common lines of sight demonstrates that 
a significant and perhaps dominant fraction of low-redshift [z < 1) Lya- forest absorption systems 
arise in extended gaseous envelopes of galaxies (Lanzetta et al. 1995a; Chen et al. 1998). But it 
has not yet been possible to extend this comparison to higher redshifts, because normal surveys 
fail to identify galaxies at redshifts much beyond z = 1. Nevertheless, determining the origin of 
high-redshift Lya-forest absorption systems bears crucially on all efforts to apply the absorbers 
as a means to study galaxies in the very distant universe and to study how the extended gas of 
galaxies evolves with redshift. Previous arguments based on the apparent low metal content and 
lack of clustering of the absorbers suggested an intergalactic origin (Sargent et al. 1980), but these 
arguments are weakened by recent Keck measurements of the actual metal content and clustering 
of the absorbers (Songaila & Cowie et al. 1996; Fernandez-Soto et al. 1996). 

To extend the comparison of galaxies and absorbers to higher redshifts, it is necessary to 
systematically identify normal galaxies at redshifts z > 1. Over the past few years, it has become 
clear that high-redshift galaxies can be reliably identified by means of broad-band photometric 
techniques. For example, several groups have determined photometric redshifts of galaxies in the 
Hubble Deep Field (HDF) to redshifts as large as z « 6 (Lanzetta, Yahil, & Fernandez-Soto 
1996; Sawicki, Lin, & Yee 1997; Connolly et al. 1997; Lanzetta, Fernandez-Soto, & Yahil 1998; 
Fernandez-Soto, Lanzetta, & Yahil 1999). Spectroscopic redshifts of nearly 120 of these galaxies 
obtained using the Keck telescope (see Fernandez-Soto et al. 1999 for a complete list of references) 
have demonstrated that the photometric redshifts are both reliable enough and accurate enough to 
establish galaxy surface density versus redshift to large redshifts (Lanzetta et al. 1998; Fernandez- 
Soto et al. 1999). 

Here we combine an empirical measure of the galaxy surface density versus redshift, obtained 
from the HDF, with an empirical measure of the gaseous extent of galaxies, obtained by Chen et 
al. (1998), to predict the number density of Lya absorption systems that originate in extended 
gaseous envelopes of galaxies versus redshift. Previous comparison of galaxies and Lya absorption 
systems at redshifts z < 1 showed that (1) galaxies of all morphological types possess extended 
gaseous envelopes and (2) the gaseous extent of galaxies scales with galaxy i?-band luminosity but 
does not depend sensitively on redshift (Chen et al. 1998). If these results apply to galaxies at 
all redshifts, then known galaxies of known gaseous extent must produce some fraction of Lya 
absorption systems at all redshifts. 

On the basis of our analysis, we find that this fraction is significant. Specifically, considering 
Lya absorption systems of absorption equivalent width W ^ 0.32 A, we find that galaxies can 
account for nearly all observed Lya absorption systems at z < 2 and that galaxies of luminosity 
Lb <; 0.05LB, can account for approximately 50% of the observed Lya absorption systems at 
higher redshifts. We further argue that if the gaseous extent of galaxies does not decrease with 
increasing redshift, then known galaxies must produce at least as many Lya absorption systems 
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as our predictions. We show that we can already explain 70% and perhaps all of the observed 
Lya absorption systems at z > 2.0 after correcting for faint galaxies that are below the detection 
threshold of the HDF images. 

We adopt a Hubble constant, Hq = 100 h km s _1 Mpc -1 , and a deceleration parameter, 
qo = 0.5, throughout the paper, unless otherwise noted. 

2. METHODS 

To calculate the predicted number density of Lya absorption systems that originate in extended 
gaseous envelopes of galaxies versus redshift, we multiply the galaxy number density per unit 
comoving volume by the absorption gas cross section of galaxies, which is known to scale with 
properties of galaxies. Our previous comparison of galaxies and Lya absorption systems at z < 1 
showed that the amount of gas encountered along a line of sight depends on galaxy impact parameter 
p and .B-band luminosity Lb but does not depend strongly on galaxy average surface brightness 
(fi), disk-to-bulge ratio D/B, or redshift z (Chen et al. 1998). Namely, galaxies of all morphological 
types possess extended gaseous envelopes, and the absorption gas cross section scales with galaxy 
.B-band luminosity but does not evolve significantly with redshift. The left panel of Figure 1 shows 
a statistically significant correlation between galaxy S-band magnitude and the residuals of the 
W versus p anti-correlation. In contrast, the right panel of Figure 1 shows no correlation at all 
between galaxy redshift and the residuals of the W versus p anti-correlation accounting for galaxy 
B-band luminosity. 

The predicted number density of Lya absorption systems n{z) originating in extended gaseous 
envelopes of galaxies can now be written 

c f°° 
n{z) = — (1 + z)(l + 2q z)- 1 ' 2 / dL B 3>(L B , z)ttR 2 (L b ), (1) 
^0 JL Bin . n 

where c is the speed of light, $>(Lb, z) is the galaxy luminosity function, and R is the gaseous extent. 
The minimum galaxy B-band luminosity i_B min is determined by the detection threshold of either 
the galaxy redshift survey or the Lya absorbing galaxy survey. It is taken to be the brighter of the 
magnitude limits at which the galaxy luminosity function or the scaling relation was determined. 
Supplementing results of Chen et al. (1998) with new measurements (Chen et al., in preparation), 
we find that galaxy gaseous radius R scales with galaxy B-band luminosity Lb as 

R_ = (LbV 
R* \L B J ' 

with 

t = 0.39 ± 0.09, (3) 
and with R* = 172 ± 27 hr l kpc at W = 0.32 A for galaxies of luminosity L B > 0.03L B , . 



(2) 
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Given a complete galaxy sample, we can evaluate equation (1) by establishing an empirical 
galaxy luminosity function without adopting any particular functional form. First, we write the 
galaxy luminosity function as a discrete sum of 5 functions, 

* (L *' «> - (£)" (1 + z) " il + 2 * 2)1/2 atthw P (Lb - Lb - ] - (4) 

where Q is the angular survey area, Da is the angular diameter distance, and L Bl is the £>-band 
luminosity of galaxy i. Next, we substitute equation (4) into equation (1) and rearrange terms to 
yield 

n(z) -l 1 T *R\L Bi ) 

where (z±,Z2) marks the boundary of each redshift bin. Equation (5) indicates that the number 
density of Lya absorption systems originating in the extended gaseous envelopes of galaxies is 
equal to the product of the galaxy surface density and the average of luminosity weighted gas cross 
sections. 



3. GALAXY SAMPLE 

Here we adopt a galaxy sample from the Hubble Deep Field (HDF) galaxy catalog published 
by Fernandez-Soto et al. (1999), which contains coordinates, optical and near- infrared photometry, 
and photometric redshift measurements of 1067 galaxies. This galaxy catalog is unique for studying 
galaxy statistics at high redshifts, because it is sensitive to objects of surface brightness down to 
ps 26.1 mag arcsec -2 and is complete to AB(8140) = 28.0 within the central 3.93 arcmin 2 area 
of the HDF image (Zone 1 of Fernandez-Soto et al. 1999). Photometric redshifts ranging through 
z ps 6 were measured from broad-band spectral energy distributions established from the optical 
(Williams et al. 1996) and infrared (Dickinson et al. 1999, in preparation) images. Spectroscopic 
redshifts of nearly 120 of these galaxies have been obtained using the Keck telescope (see Fernandez- 
Soto et al. 1999 for a complete list of references). Comparison of the photometric and spectroscopic 
redshifts shows that the photometric redshifts are accurate to within an RMS relative uncertainty 
of Az/(l + z) < 0.1 at all redshifts z < 6 (Lanzetta et al. 1998; Fernandez-Soto et al. 1999), except 
that there are not spectroscopic redshifts available for comparison at redshifts z = 1.5 to z = 2.0. 
Reliable galaxy statistics as a function of redshift can therefore be measured on the basis of the 
HDF galaxy catalog. 



4. ANALYSIS 

The goal of the analysis is to compare the predicted number density of Lya absorption systems 
originating in extended gaseous envelopes of galaxies (using the HDF galaxy catalog) with the 
observed number density of Lya absorption systems. 
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Here we evaluate equation (5) for galaxies observed in the HDF, given the photometric redshift 
measurements. Galaxy £?-band luminosity is calculated from the apparent F814W magnitude cor- 
rected for luminosity distance, color k correction (between rest-frame B band and observed-frame 
F814W band, calculated based on the spectral type determined from the photometric redshift tech- 
niques), and bandpass k correction. Errors in the predicted number density are estimated using a 
standard bootstrap method, which incorporates the sampling error with the uncertainties of pho- 
tometric redshift measurements and the scaling relation. To model the uncertainty of photometric 
redshift measurements correctly, Lanzetta et al. (1998) have shown that both photometric error 
and cosmic variance with respect to the spectral templates must be accounted for. The first can 
be simulated by perturbing galaxy photometry of different bands within the photometric errors. 
The second can be characterized by the RMS dispersion of photometric and spectroscopic redshift 
measurements, which is 0.08 at z < 2 and 0.32 at z > 2 (Lanzetta et al. 1998). The uncertainty of 
photometric redshift measurements is then calculated by forming a quadratic sum of the two. 

The results are shown in Table 1 and Figure 2 in comparison with observations. In Figure 2, 
filled circles represent the observations and open circles represent our predictions. The observations 
at z < 1.5 are taken from the sample 2 of Weymann et al. (1998), which we normalize to refer to Lya 
absorption lines of W ;> 0.32 A. The observations at z > 1.6 are taken from the sample 9 of Bechtold 
(1994)0. We select a bin size, Az = 0.5, when calculating the predicted values. Experiments with 
different bin sizes show that the redshift distribution of the predicted number density of Lya 
absorption systems is insensitive to the selected bin size. We also repeat the calculation, using a 
deceleration parameter, qo = 0.0. It turns out that the predicted number density of Lya absorption 
systems does not depend sensitively on the adopted deceleration parameter either. 

Our predictions have so far been limited to including galaxies of luminosity Lb > Ls min - 
To estimate the contribution of faint galaxies to the predicted gas cross section, we calculate the 
number density of Lya absorption systems originating in galaxies of luminosity Lb < Ls min by 
choosing a Schechter luminosity function. We evaluate equation (1) by adopting a faint-end slope 
obtained by Ellis et al. (1996), which is a = — 1.41^oo7i an d by extrapolating the scaling relation 
(equation 2 and 3) to obtain the extended gaseous radii for galaxies of luminosity ^ Lb < 0.03Lb, ■ 
Increments of the predicted number density versus redshift due to the inclusion of faint galaxies 
are shown in Figure 2 as well (crosses with dashed horizontal bars). 



3 Note that the data from Weymann et al. included all the Lya absorption lines with or without associated metal 
absorption lines, while the data from Bechtold included only the Lya absorption lines without associated metal 
absorption lines. However, the correction for including metal-line associated systems in the high-redshift sample is 
expected to be approximately a few percent (Bechtold 1994; Frye et al. 1993). In comparison to the large uncertainty 
of the number density of high-redshift Lya absorption systems, it is clear that the absence of metal-line associated 
systems in the high-redshift sample does not change the redshift distribution of Lya absorption systems by a noticeable 
amount. 
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5. DISCUSSION 

The HDF images have provided us with a unique chance to relate statistical properties of high- 
redshift galaxies to statistical properties of Lya absorption systems, thereby studying the origin 
of Lya absorption systems at high redshifts. Combining an empirical measure of galaxy surface 
density versus redshift with an empirical measure of gaseous extent of galaxies, we have predicted 
the number density of Lya absorption systems that originate in extended gaseous envelopes of 
galaxies versus redshift. 

For Lya absorption systems of absorption equivalent width W ;> 0.32 A, comparison of the 
predicted and observed number densities of Lya absorption systems shows that (1) known galaxies 
with known gas cross sections can account for all the observed Lya absorption systems to within 
measurement errors at redshifts < z < 1, (2) known galaxies of luminosity Lb ^ 0.03 with 
unevolved gas cross sections can account for all Lya absorption systems at redshifts 1 < z < 2, 
and (3) known galaxies of luminosity Lb ^ 0.05 Lb, with unevolved gas cross sections can account 
for approximately 50% of observed Lya absorption systems at z > 2. Apparently, bright galaxies 
alone can account for a significant portion of the observed Lya absorption systems at all redshifts. 
After correcting for faint galaxies that are below the detection threshold of the HDF image, we can 
already explain 70% and perhaps all of the observed Lya absorption systems at z > 2.0 to within 
errors. Here we discuss three factors that may invalidate the result. 

First, we consider evolution of the galaxy luminosity function. There has been some evidence 
showing that galaxies at the faint end of galaxy luminosity function evolves significantly with 
redshift (e.g. Ellis 1997, and references therein). Namely, faint galaxies may be more numerous at 
higher redshifts. To address this issue, we evaluate equation (1) with ^B min ~ to calculate the 
relative change in the predicted number density, given a different faint-end slope. It turns out that 
the predicted number density of Lya absorption systems increases by as much as 0.3 dex, if we 
vary the faint-end slope from a = —1.4 to a = —1.6. Apparently, a steeper faint-end slope brings 
the predicted values closer to the observed ones. We conclude that known galaxies must produce 
at least as many Lya absorption systems as the predictions. 

Next, we consider evolution of neutral gas surrounding galaxies. In the right panel of Figure 
1, we show that the residuals of the W versus p anti-correlation after accounting for the scaling of 
galaxy Z?-band luminosity do not correlate with galaxy redshift. Although there is no direct 
observational evidence to support the extension of the redshift independence of the scaling relation 
to redshifts beyond z = 1, the gaseous extent of galaxies is unlikely to decrease with increasing 
redshift in these early epochs. Theoretically, it is believed that galaxies were formed through 
accretion of cooled gas over time, indicating a larger gaseous extent of galaxies at higher redshifts. 
Observationally, the mass density of neutral gas measured from damped Lya absorption systems 
increases significantly with increasing redshift at redshifts between z = 1.6 and z = 3.5 (Lanzetta, 
Wolfe, & Turnshek 1995; Storrie-Lombardi, McMahon, & Irwin 1996), supporting that galaxies may 
possess more neutral gas and therefore no less gaseous extent at higher redshifts. We again conclude 
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that known galaxies must produce at least as many Lya absorption systems as the predictions. 

Finally, we consider galaxy clustering effects. Due to limited resolutions of spectroscopic 
observations, Lya absorption lines with small velocity separation are sometimes blended together 
and considered as one absorption system. Lanzetta, Webb, & Barcons (1996) first reported an 
absorption system which may arise in a group or cluster of galaxies. It was later confirmed by Ortiz- 
Gil et al. (1999), who identified eight absorption features in this system based on a spectrum of 
higher spectral resolution. Their analysis suggested that the degree of clustering of Lya absorption 
systems may be underestimated on velocity scales of several hundred kilometers per second (see also 
Fernandez-Soto et al. 1996). Similarly, the number density of Lya absorption lines derived from 
a known galaxy population would be in excess of the observed number density of Lya absorption 
lines, because galaxies are strongly clustered and the predicted number density does not suffer from 
the line blending effect. As a result, the fraction of Lya absorption systems originating in extended 
gaseous envelopes of galaxies may be overestimated. 

To estimate the amount of excess at different redshifts, we calculate the expected number of 
neighbouring galaxies of luminosity Lb > ^B min , -^neighbours, that are seen at impact parameters 
p <z 200 kpc along a line of sight with a velocity difference of Av from an absorption redshift. We 
adopt a two-point correlation function measured by Magliocchetti & Maddox (1999) and Arnouts 
et al. (1999) for the HDF galaxies at redshifts <^ z <^ 4.8, which is characterized by 

£(r,z) = ( r /r )^(l + z)-( 3+£ ) (6) 

with 7 = 3 + e = 1.8 and r$ ~ 1.7 h^ 1 Mpc, and evaluate the integral of the two-point correlation 
function over a comoving volume spanned by 5v in depth and 200 kpc in radius. The velocity 
span, Av, is taken to be the spectral resolution of which the observations were carried out. It is 
250 km s _1 for the sample from Weymann et al. (i.e. absorption systems at z < 1.5) and 75 km s _1 
for the sample from Bechtold (i.e. absorption systems at z > 1.5). We show in Table 2 that the 
estimated excess in our predictions (represented by the number of neighbouring galaxies) can be at 
most a factor of four at z < 2, but is negligible at z > 2. Therefore, we conclude that the fraction 
of Lya absorption systems originating in extended gaseous envelopes of galaxies remains significant 
even after correcting for galaxy clustering. 

In summary, we present our first attempt at relating statistical properties of galaxies to sta- 
tistical properties of Lya absorption systems at high redshifts based on the HDF galaxy catalog. 
Combining an empirical measure of galaxy surface density versus redshift with an empirical measure 
of gaseous extent of galaxies, we have predicted the number density of Lya absorption systems that 
originate in extended gaseous envelopes of galaxies versus redshift. We show that at approximately 
50% and as much as 100% of observed Lya absorption systems of W ;> 0.32 A can be explained by 
extended gaseous envelops of galaxies. The result remains valid after taking into account possibile 
evolutions of faint galaxies and extended gas around galaxies, as well as galaxy clustering effects. 
Therefore, we conclude that known galaxies of known gaseous extent must produce a significant 
fraction and perhaps all of Lya absorption systems over a large redshift range. 
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Table 1. NUMBER DENSITY OF Lya ABSORPTION SYSTEMS WITH W £ 0.32 

A— PREDICTED AND OBSERVED 



n(z) 



Redshift 


Predicted 


Observed 


Lmin/ Lb* 


0.0 < z < 0.5 


13 


27 


0.03 


0.5 < z < 1.0 


29 


28 


0.03 


1.0 < z < 1.5 


27 


30 


0.03 


1.5 < z < 2.0 


34 


30 


0.03 


2.0 < z < 2.5 


26 


58 


0.05 


2.5 < z < 3.0 


25 


70 


0.09 


3.0 < z < 3.5 


24 


77 


0.12 


3.5 < z < 4.0 


13 


89 


0.17 



a Here we adopt an Lb* corresponding to Mb* = 
-19.5 (Ellis et al. 1996). 



Table 2. ESTIMATED EXCESS OF PREDICTED Lya ABSORPTION SYSTEMS 



Redshift 


-^neighbours 


0.05 < z < 0.5 


1.7 


0.5 < z < 1.0 


2.7 


1.0 < z < 1.5 


3.3 


1.5 < z < 2.0 


1.1 


2.0 < z < 2.5 


0.8 


2.5 < z < 3.0 


0.6 


3.0 < z < 3.5 


0.5 


3.5 < z < 4.0 


0.4 
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Fig. 1. — Left: The residuals of the W vs. p anti-correlation as a function of galaxy i?-band 
magnitude. Right: The residuals of the W vs. p anti-correlation accounting for galaxy f?-band 
luminosity as a function of galaxy redshift. Solid circles represent elliptical or SO galaxies, triangles 
represent early-type disk galaxies, and squares represent late-type disk galaxies. Points with arrows 
indicate no detections at a 3-a upper limit. The coefficients, a and (3, were determined from a 
maximum likelihood analysis previously discussed in Chen et al. (1998). 



-12- 



~~\ 1 r 



"i 1 1 1 1 1 1 1 1 r 



W > 0.32 A 



N 



O 



1.5 



i — & 



& — i 



— 4 



« 1 



I I I I I I I I I I I I I I L 

0.2 0.4 0.6 

log (1 + z) 



0.8 



Fig. 2. — Comparison of the observed number densities of Lya absorption systems of absorption 
equivalent width W > 0.32 A (closed circles) and the predicted ones produced in extended gaseous 
envelopes of galaxies observed in HDF (open circles). The absorption line observations are from 
Bechtold (1994) for data at z > 1.6 and from Weymann et al. (1998b) for data at z < 1.5. We have 
multiplied the data from Weymann et al. by a factor of exp ( ' Q ^ 73 — J in order to be compared with 
observations from Bechtold. The vertical bars represent 1-a errors in the observed and predicted 
values and the horizontal bars indicate the adopted bin size, Az = 0.5. Crosses indicate the 
predicted number densities of Lya aborption systems arising in galaxy gaseous envelopes after 
correcting for galaxies fainter than the detection limit of the HDF image by taking Ls min = 0.0 
and adopting a galaxy luminosity function obtained by Ellis et al. at 0.15 < z < 0.35 (1996). 



